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We use an effective QCD model (ladder-QCD) to explore the phase diagram for 
chiral symmetry breaking and restoration at finite temperature with different u, d 
quark chemical potentials. In agreement with a recent investigation based on the 
Nambu-Jona-Lasinio model, we find that a finite pion condensate shows up for high 
enough isospin chemical potential fij = (fi u — F° r small fij the phase diagram 

in the , T) plane shows two first order transition lines and two critical ending 
points. 

PACS numbers: ll.10.Wx, 12.38.-t, 25.75.Nq 



I. INTRODUCTION 

In the last few years the study of QCD at finite density has become rather important. 
In particular it has been established that at zero temperature and in the high density limit 
a color superconducting phase exists (for a review see ^|). From a phenomeno logical point 
of view there are two areas where finite density is relevant. The first area is the realm 
of compact stellar objects where the central density can reach values up to ten times the 
saturation density p, with p « 0.14 fm -3 evaluated as the inverse of the volume of a sphere of 
radius 1.2 fm. Since the temperature of a compact star is much smaller than the typical color 
superconducting gap (of order of tens of MeV) one can safely consider the limit T — > 0. The 
second area can be found in heavy ion physics. However, in this case, color superconductivity 
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is not relevant given the large entropy per baryon produced in heavy ion collisions. But 
lere another important feature of QCD might be relevant. According to different models 
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the phase diagram of QCD in the plane (/xb,T) exhibits a tricritical point. Since 
this point should be located at moderate density and temperature there is some possibility 
of observation in heavy ion experiments. Another important point in heavy ion physics is 
the fact that in the experimental setting there is a non zero isospin chemical potential, fij. 



Studies at finite fij have been the object of several papers ((J Q, His, 10, 11] but mainly in 
the regime of low temperature and high baryon chemical potential. The first complete study 
of the phase diagram in the three-parameter space (/ib, T, fii) has been made in the context 
of a Random Matrix model [l2|. It has been found that the first order transition line ending 
at the tricritical point of the case /if = actually splits in two first order transition lines and 
correspondingly two crossover regions are present at low values of baryon chemical potential. 
The existence of this splitting has also been shown in the context of a Nambu-Jona-Lasinio 
m0 del in Q. It S hon,d also be noticed that in Q the NJL m ode, ha S been a Ugm ented h y 
the four-fermi instanton interaction relevant in the case of two flavors. These authors have 
found that the coupling induced by the instanton interaction between the two flavors might 
wash completely the splitting of the first order transition line. This happens for values of 
the ratio of the instanton coupling to the NJL coupling of order 0.1-0.15. 

In this paper we will consider the effect of a finite isospin chemical potential in a model 
(ladder-QCD) where the existence of a tricritical point was shown several years ago [2]. The 
reason of doing this analysis in a model different from the NJL model is due to the fact that 
QCD at finite baryon density is difficult to be studied on the lattice (however it should be 
noticed that recently a new technique has been proposed 3] and a first evaluation of the 
tricritical point has been given in |l6j). It is therefore important to study certain features in 
different models in order to have a feeling about their universality. For instance the existence 
of a tricritical point seems to enjoy such a characteristic. What we are presenting here is 
a preliminary study, and therefore we will restrict the analysis at small isospin chemical 



potential. The interest for this topic is due to results from lattice simu 
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ations and effective 



18, 



19, 201. We also 



theories which show the existence of a phase transition at finite fij 
ignore the effects from color superconductivity, since these are present only for temperatures 
lower than some tens of MeV. 

In our model all the three flavors are present, and the relevant instanton effects would 
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give rise to a six-fermi contact interaction. Therefore it is not clear if these effects will wash 
out the splitting as in the two flavor case 3]. We will consider this problem in a future 
work. 



II. THE MODEL (REVISITED) 



In this Section we will review a model that was used several years ago to describe the 



chiral phase of QCD both at zero temperature |21|, |22| and at finite temperature and density 
j3, 13]. This model is an approximation to QCD based on the evaluation of the effective 
potential at two-loop level and on a parametrization of the self-energy consistent with the 
OPE results. The effective action that we evaluate is a slight modification (see for instance 
P]) of the Cornwall- Jackiw- To mboulis action for composite operators I^Zl^. We recall here 
the major steps of this calculation. We start from the Cornwall- Jackiw- Tomboulis formula 



r[s] 



T 2 [S] + Tr 



s 



.5T2 
5S 



Tr In 



s ' + 



6T2 
5S 



where the free fermion inverse propagator is 



+ counterterms 



S 1 (p) = ip - m (2) 

and T 2 is the sum of all 2PI diagrams with propagator S, which has to coincide with the 
exact fermion propagator at the absolute minimum of T. Thus 5* is the dynamical variable 
in this variational approach. However it turns out useful to trade S for 

^r 2 



5S 



(3) 



which coincides with the fermion self-energy at the minimum of T. 

In the present model (ladder-QCD) we will make the very rough approximation of evalu- 
ating r 2 at the lowest order. That is, we evaluate r 2 at two-loops with one gluon exchange. 
The relevant Feynman diagram is given in Fig. ^ It turns out that this approximation 
works rather well phenomenologically (see for instance 2^). Therefore, at this order the 



effective action is simply 



22 



r[E] = - Tr In (Sq 1 - S) - - Tr (S E) + c.t. 
= - Tr In (Sq 1 - S) + r 2 [E] + c.t. 



(4) 
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FIG. 1: The two-loop diagram needed to evaluate T%. 



As in refs. 2l|,|22j], we use the following parametrization for S 



S{p) = iA{p)p + B(p) + i7 5 C(p) (5) 

Then, by working in the Landau gauge, it is possible to show that no renormalization of 
the wave function is required and also that the Ward identity at this order is satisfied by 
taking the free quark-gluon vertex and the free gluon propagator. We also consider the 
so-called rigid case, where the strong coupling g is considered fixed at p 2 = M 2 , where M 
is a convenient mass scale to be fixed later. In this way the relation between the scalar and 
pseudoscalar contribution to S and the terms B, C in Eq. (j5J) can be easily inverted, leading 
to the following expression for the effective action, fully expressed in terms of E (|22|) 

r[s] = r 2 p] + r log [s] (6) 

where, separating S = S s + ij^ p 

T 2 [S] = ~ 8 -^f J ^tr [E s (g 2 )D S s (g 2 ) + S p (g 2 )D S p (g 2 )] (7) 

with f2 4 being the four- volume and C2 = 4/3 for N c = 3 the quadratic Casimir of SU(3) C . 
Besides M, also g is a parameter of the model. The one-loop term is 

r log [E] = -Tr log (So" 1 - E) = -Tr log [ip-m- £ s (p 2 ) - i^ P (p 2 )] (8) 

where the scalar and pseudoscalar parts of the dynamical variable are matrices in SU(3) 
flavor space (as well as m), related to the scalar and pseudoscalar quark condensates through 
the following equation 



S s (p 2 ) + i7 5 S p (p 2 ) = + ij 5 p) f(p 2 ) 



(9) 
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The function f(p 2 ) which contains the momentum dependence of the self-energy will be 
discussed in a moment. The fields 



(■Sab) 



(Pab) 



3C 2 g 2 
ANM 3 

SC 2 g 2 



(^a^b) 



(* a i75*6> 



(10) 



(11) 



ANM 3 

will be determined by minimizing the effective action. In Eq. (JSJ), m is the mass matrix in 
flavor space which is taken diagonal 

n n \ 



m 



\ 



m u 











m d 











m s 



J 



The function f(p 2 ) will be chosen requiring that it goes to a constant for p — > and as 1/p 2 
(mod log terms) for large values of p as suggested by the OPE expansion. By introducing a 
dimensionless variable x 2 = p 2 /M 2 we will consider the following family of functions 



fN{x 2 

M 



1 + X + X + ... + X 



1 + X 2 + X 4 + 



r 2N-2 



In the limit N — > oo we get the function used in 

f(x 2 ) 



21 



22 



+ x 



2N 



M 



Notice that for x — > we get 



-_e(i-x 2 ) + e(x 2 

f N {x 2 ) 



M 



2iV 



(12) 



(13) 



(14) 



Now, let us consider for simplicity the chiral limit and zero chemical potentials. In this case 
it is simple to get the mass-shell condition from the one loop term in Eq. (jHJ) (see for instance 

B) 

P 2 + (s) 2 f N = (15) 

where (s) is the field proportional to the scalar condensate (see Eq. (fTO])). If we want to 
recover, at least in the infrared regime, a free particle-like dispersion relation (as for instance 
happens in four fermion theories), we see from Eqs. (fUjl and (fTKj) that we need N > 2. In 
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this paper we will choose N = 2. Notice that the choice N = 1 would lead to the following 
dispersion relation in the limit of small momenta 

p 2 (l-2(s) 2 ) + M 2 (s} 2 + ... = (16) 

This might give rise to problems in the broken phase where the coefficient of p 2 could 
become negative. However no difficulties arise for the determination of the critical points 
where (s) ~ 0. On the contrary the equation of state could be affected. 
We can thus evaluate explicitly the effective potential 

v - k < 17 > 

which is UV- finite in the chiral limit, whereas it needs to be properly renormalized in the 
massive case . We have employed the following normalization condition 

dV 

d(m a (y a y a }) min ~ 1 ^ 

In the chiral limit, this requirement is equivalent to the Adler-Dashen relation (see for 
instance [s^l). Here we will require the validity of this equation at the values of the quark 
current masses. 
By defining a a = m a /M and 

Xa = -^(Va) (19) 

the normalization condition, using Eq. (jl(Jj) . can be written as 

dV a 3M 4 _ 3M 4 

d(^a) m ln ~~^ C ~~J{M) (20) 

where we have introduced the parameter 

(21) 

The chemical potential and the temperature dependence are introduced following standard 
methods 



26| (see for example |2j). In particular the chemical potential is introduced from 
the very beginning via the usual substitution p v —>■ (po + ijti, p) in p appearing in the Dirac 
operator in Eq. (jSJ). On the other hand the temperature dependence is introduced by 
substituting to p the Matsubara frequency u n = (2n + l)nT in all the p dependent terms 
appearing in the effective action. The reason for this asymmetrical treatment is that the p 
dependence in the self-energies becomes relevant only at p ^> M whereas, as we shall see we 
will be interested in chemical potentials lower than M. 
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III. RESULTS AT FINITE TEMPERATURE AND DENSITY 

In order to get the effective action we need to calculate the determinant of the operator 
appearing in Eq. (JEJ). We set to zero the strange quark chemical potential p s = and define 
Pi = (Pu — Pd)/2 and \ib = (p u + Pd)/2- The operator is given by the following 3 ® 3 matrix 
in flavor space: 



'i(u n + ip u )lQ + ip '■ 7 - K p/275 ^ 

-p/275 i(w n + ipdho + ip ■ 7 - F d (22) 

y iu^ + ip ■ j -F s j 

where we have defined F a = m a + f2Xa] a = u,d,s and with f 2 given in Eq. (|T2j) with 
N = 2. Also p 2 = + \p\ 2 . Here p is related to the charged pion condensate 

P = ~Sp ({Ul5d) ~ (J75M)) (23) 

We have directly set to zero the hyper-charged condensates in the strange sector, since 
we do not expect the formation of such condensates for /i s = 0. The strange sector thus 
factorizes and the determination of the phase diagram for approximate chiral symmetry 
restoration is performed by studying the behavior of the light (uu) and (dd) quark conden- 
sates, independently on the strange quark condensate. To evaluate the effective action, we 
perform the sum over the Matsubara frequencies which solve the mass-shell condition given 
by the vanishing of the determinant in Eq. (J22j) using standard methods Although 
formally straightforward, the calculation is a very hard numerical task. Actually, at each 
integration step on \p\, we have to solve a twentieth order algebraic equation in u n in the 
u — d sector. The part relative to the strange quark is obviously easier to deal with. The 
coefficients of this equation depend also on the parameters Xu, Xd, p, Pu, Pd, m u , fna- 

Before discussing the results at finite density and temperature let us review how we fix 
the parameters appearing in the effective action. This is done by looking at T = [l = 0. The 
parameters that we have to fit are c, M, m u , ma, m s . These are obtained by using as input 
parameters the following physical quantities m n ±, m#-±, m^o, f n and fx- The results of the 
fit are given in Table I, whereas the values of the input experimental quantities, together 
with the result we get from the fit procedure are given in Table II. 
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Parameters 


Fitted Values 


M 


529 MeV 


c 


1.0 


m u 


4.4 MeV 


m d 


6.2 MeV 


m s 


110 MeV 



TABLE I: Fit of the parameters. 



Input parameters 


Fitted Values 


Experimental Values 


m^± 


139 MeV 


139.6 MeV 


m K ± 


494 MeV 


493.7 MeV 


m K o 


499 MeV 


497.7 MeV 


u 


92 MeV 


92.4 MeV 


fK 


105 MeV 


113 MeV 



TABLE II: Comparison between the values of the input parameters as obtained from the fit and the 
experimental results. 

With these values of the parameters, we find that at T = // = the quark condensate in 
chiral limit has the value 

(^) = -(248 MeV) 3 (24) 

whereas in the massive case 

(mi) = -(251 MeV) 3 , (dd) = -(253 MeV) 3 , (ss) = -(305 MeV) 3 (25) 
and, by defining the constituent quark masses a la Politzer j^] 

M const = £(p 2 = AMj onst ) (26) 

we get 

M s = 385 MeV, M u4 = 256 MeV (27) 

where, here and in the following the light quarks have been taken degenerate with mass 
rh= (m u + m d )/2 = 5.3 MeV. 
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From the general study at T = fi = and in the chiral limit, one finds also that in order to 
break the chiral symmetry one must have 

c < 1.37 (28) 

This condition is satisfied by our choice of parameters. 

When ix j = it was shown that in the model discussed in [2 1 ] there is a tricritical point in 
the chiral limit. The model we are presenting here is essentially the same model with some 
slight modifications, as for instance, the choice of the function f'2{x). However the tricritical 
point is still present as it can be seen from Fig. El (central line), obtained with the choice of 
parameters of Table I. 



T(MeV) 




50 100 150 200 250 300 



350 
|Ll B (MeV) 



FIG. 2: Phase diagram for chiral symmetry in the (hb,T) plane. For fij = (central line), the 
cross-over transition line starts from the point Pt = (0, 135) and ends at the point E = (240,97). 
The line between E and the point P^ = (325, 0) is the line for the first order transition with 
discontinuities in the (uu) and {dd) condensates. For llj = 30 MeV (side lines), the two cross- 
over transition lines start from the point Pj- = (0, 135) and end at the points E u = (210, 97) and 
Ed = (270,97). The lines between E u and the point Pu = (295,0) and between Ed and the point 
P^ = (355, 0) are the lines for the first order transitions with discontinuities in the {uu) and (dd) 
condensates respectively. 

A complete analysis of the full three parameter space (lib,T, hi) has not yet been com- 
pleted especially in relation with the pion (and hyper-charged) condensate. We have exam- 
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ined the case T = fig = and we have found that there is a phase transition indicated by 
a finite pion condensate starting at fij = 70 MeV. Since in our model m n ~ 140 MeV we 



We found also that for 



agree with the results found in the literature 
m = the critical fii is zero. Therefore we will limit our considerations at small isospin 
chemical potential, say pi = 30 MeV, where we expect the pion condensate p to vanish. In 
this situation the determinant in Eq. (}2*2*|) factorizes and the effective action is given by a 
sum of three independent terms, one for each flavor. Therefore the action is the same as for 
Hi = 0, with each flavor evaluated at its own chemical potential. It follows that the light 
flavor terms show the same tricritical structure exhibited from the central line in Fig. El 
for fii = 0. Consequently the phase diagram we obtain for a small, fixed isospin chemical 
potential (/!/ = 30 MeV) is described by the two side lines in Fig. El with each flavor u 
and d showing the same structure as the central line but with a split chemical potential 
Hu,d = i Hi. Notice also that although the Figure extends up to zero temperature, this 
part should not be taken too seriously due to the existence of color superconductivity. 

IV. CONCLUSIONS 

In this paper we have discussed an approximate model of QCD (ladder-QCD) at finite 
temperature and densities. In particular we have considered the experimentally important 
situation of a non-vanishing isospin chemical potential. This situation has been already 
explored 



121 and 



jy various authors previously and we confirm, in particular, the results found in 
13j about the split ting of the first order transition line in the plane (hb,T), for 
small Hi. As pointed out in [l3j] this result could be relevant for ion physics since the first 
order transition line is split symmetrically with respect to the original line at \ii = 0. This 
implies a reduction of the value of the baryon chemical potential at the tricritical point of an 
amount given by /i/, making easier the possibility of discovering it experimentally. Since it 
is very difficult to perform first principle analysis of QCD at finite baryon density, we think 
that it is important to show that certain features as the existence of the tricritical point and 
the possible splitting of the first order transition line are common to several models. This 
suggests that these features might possess some universal character. However we should 
stress that the result of the splitting of the first order transition line is strictly related to the 
factorization in flavor space. For instance, in the two flavor case, the four-fermi interaction 
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due to the instanton effects leads to a mixing of the flavors that, if sufficiently large, might 
wash out the mixing. We think that this point needs further analysis in the more complete 
scheme with three flavors. 
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